Generalized synchronization: a modified system approach* 



Alexander E. Hramo\0 and Alexey A. Koronovskif] 

Faculty of Nonlinear Processes, Saratov State University, Astrakhanskaya, 83, Saratov, ^10012, Russia 

(Dated: February 8, 2008) 

The universal mechanism resulting in the generalized synchronization regime arising in the chaotic 
oscillators with the dissipative coupling has been described. The reasons of the generalized synchro- 
nization occurrence may be clarified by means of a modified system approach. The main results 
are illustrated by unidirectionally coupled Rossler systems, Rossler and Lorenz systems and logistic 
maps. 
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Chaotic synchronization is one of the fundamental phe- 
nomena, widely studied recently 0, having both the- 
oretical and applied significance (e.g., for information 
transmission by means of deterministic chaotic signals 
0,0], in biological and physiological Q tasks, etc.). Sev- 
eral different types of chaotic synchronization of coupled 
oscillators, i.e. generalized synchronization (GS) [a, III, 
hose synchronization (PS) y|, lag synchronization (LS) 
and complete synchronization (CS) @ are well known. 
There are also attempts to find unifying framework for 
chaoti c sy nchronization of coupled dynamical systems 

One of the interesting and intricate types of the syn- 
chronous behavior of unidirectionally coupled chaotic os- 
cillators is the generalized synchronization. The pres- 
ence of GS between the response x r (i) and drive Xd(t) 
chaotic systems means that there is some functional re- 
lation x r (t) = F[xd(t)] between system states after the 
transient finished. This functional relation F[-] may be 
smooth or fractal. According to the properties of this 
relation, GS may be divided into the strong synchroniza- 
tion and week synchronization, respectively |l3j |. There 
are several methods to detect the presence of GS be- 
tween chaotic oscillators, such as the auxiliary system 
approach y~4| or the method of nearest neighbors [a, 0] • 
It is also possible to calculate the conditional Lyapunov 
exponents (CLEs) to detect GS. The regimes of 

LS and CS are also the particular cases of GS. 

This paper aims to explain GS arising. We show the 
physical reasons leading to GS appearance in unidirec- 
tionally coupled chaotic systems. The causes of the 
generalized synchronization arising may be clarified by 
means of a modified system approach. 

Let us consider the behavior of two unidirectionally 
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coupled chaotic oscillators 

-k d = H(x d ,g d ) 

x r = G(x r ,g r ) + eP(x rf ,x r ), 



(1) 



where Xd, r are the state vectors of the drive and response 
systems, respectively; H and G define the vector fields 
of these systems, g^ and g r are the controlling param- 
eter vectors, P denotes the coupling term and e is the 
scalar coupling parameter. If the dimensions of the drive 
and response systems are Nd and N r respectively, the 
behavior of the unidirectionally coupled oscillators (JIJ is 
characterized by the Lyapunov exponent (LEs) spectrum 
Ai > A2 > • • • > A Nd+ N r ■ Due to the independence of the 
drive system dynamics on the behavior of the response 
one, the Lyapunov exponent spectrum may be divided 
into two parts: LEs of the drive system Af > • • ■ > Afj 
and the CLEs [H E3 X[ > ■ ■ ■ > \ r N ■ The condition of 
GS is X\ < (see [l3| for detail). 

The GS manifestation is mostly considered for two 
identical systems with equal or mismatched parameters 
and diffusive type of unidirectional coupling. Therefore, 
let us consider such systems first, while the case of dif- 
ferent systems and others coupling types will be briefly 
discussed later. In the case of identical systems the di- 
mensions of the drive and response oscillators are equal 
(Nd = N r = N) and the equations may be rewritten 
as 



x d = H(x d ,g d ) 

x r = H(x r ,g r ) + eA(x d - x r ), 



(2) 



where A = {tfy} is the coupling matrix, 8a — or 1 and 
Sij = (i 7^ j). It is clear, that the dynamics of the re- 
sponse system may be considered as the non- autonomous 
dynamics of the modified system 



x?7i H (x m ,g r ,s) 
under the external force eAx^ 

x m = H'(xm, g r , e) + eAx d , 



(3) 



(4) 



where H'(x, g) = H(x, g) — eAx. Note, that the 
term —eAx brings the dissipation into the modified sys- 
tem ©. 
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So, GS arising in J5J with parameter e increasing may 
be considered as a result of two cooperative processes tak- 
ing place simultaneously. The first of them is the growth 
of the dissipation in the system © and the second one 
is the increasing of the amplitude of the external signal. 
Both processes are correlated with each other by means 
of the parameter e and can not be realized in the cou- 
pled oscillators system @ independently. Nevertheless, 
let us consider these processes separately to understand 
better the mechanisms of GS arising. We start our con- 
sidering with the autonomous dynamics of the modified 
system 

For this modified system, x m (f), the quantity e is the 
dissipation parameter. When e is equal to zero the dy- 
namics of the modified system x m (i) coincides with the 
response system x r (i) without coupling. With increas- 
ing of the dissipation parameter e the dynamics of the 
modified system (j3Jl should be simplified. Therefore, the 
system x m (i) has to undergo a transition from chaotic 
oscillations to periodic ones, and, perhaps, to the sta- 
tionary state (if the dissipation is large enough). In this 
case one of the Lyapunov exponent, A™, of the modified 
system is equal to zero (or negative if the stationary state 
takes place), all other Lyapunov exponents are negative 
(0 > A™ > ■ • ■ > A™_ x ). It is important to note, that 
the Lyapunov exponent spectrum of the modified sys- 
tem © differs from CLE spectrum A^ > • • ■ > A^y of the 
system |2"|. as CLEs depend on the drive system dynam- 
ics in contrast to Lyapunov exponents of the modified 
system x m (t). Therefore, nobody can draw a conclusion 
about the appearance of GS in the coupled oscillators 
system (J5J| taking into account only Lyapunov exponents 
of the modified system ©. 

On the contrary, the external signal in tends to im- 
pose the dynamics of the drive chaotic oscillator x d (t) on 
the modified system x m (i), and, correspondingly, com- 
plicate its dynamics. Obviously, GS may take place 
only if proper chaotic dynamics of the system x m (t) is 
suppressed by the dissipation. Only under this con- 
dition the current state x m (t) of the modified system 
will be determined completely by the external signal, i.e. 
x m (£) = F[xd(i)]. According to the equation 10} , the 
functional relation x r (f) = F[xd(i)] between the response 
and drive systems will also take place, and, therefore, GS 
will be observed. 

So, GS arising in the system (J2J is possible for such 
values of e parameter when the modified system x m (i) 
demonstrates the periodic oscillations or the stationary 
state. It is well known, that even the harmonic external 
signal can cause the chaotic oscillations in the dynamical 
system with periodical dynamics. Therefore, the periodic 
regime should be stable enough for the external force not 
to excite proper chaotic dynamics of the modified system. 
So, the difference between the parameter values s p when 
the periodic oscillations take place in the system and 
£cs when GS in the system can be observed has to 
be large enough. At the same time, the amplitude of the 
external signal is small enough in comparison with the 



amplitude of periodic oscillations in the modified system 
x m (t) (in the case when the periodical regime takes place 
in J3J|). So, the generalized synchronization looks like the 
week chaotic excitation of the periodical motion. 

The similar conclusion is also correct for the station- 
ary regime of the system x m (t) when GS manifests itself 
as chaotic perturbation of the fixed state. The system 
dynamics can be considered as transient converging to 
the "fixed" point moving under the external force in the 
phase space of the modified system ©. Let us suppose 
now that controlling parameters g r)( j of the considered re- 
sponse and drive systems differ from each other slightly 
and the value of parameter e is large enough. In this 
case the transient is very short and the state of the modi- 
fied system follows the perturbed "fixed" state essentially 
small time r of delay, therefore the regime of LS can be 
observed. 

Let us consider several examples of GS to illustrate 
the concept described above. As the first system we have 
selected two unidirectionally coupled Rossler oscillators 

id = -UdVd - Zd, 
ijd = uj d x d + ay d , 
z d =p + z d (x d -c), 

(5) 

— X r I, 

ij r = ui r x r + ay r , 
z r = p + z r (x r - c) , 

where e is a coupling parameter, u) r = 0.95. The control 
parameter values have been selected by analogy with 
as a = 0.15, p = 0.2, c = 10.0. Correspondingly, the 
modified Rossler system is 

x m oj r y m z m sx lni 

y m = u> r x m + ay m , (6) 

z m — p ~\~ z ra {x ra c) . 

In Fig I, a the bifurcation diagram for the system @ is 
shown. It is clear that this system undergoes the tran- 
sition from chaotic to periodic oscillations through the 
inverse cascade of period doubling. The dependence of 
two largest Lyapunov exponents A™ x on the parameter e 
is presented in Fig. 1,6. One can easily see, that starting 
from the value e p w 0.06 the periodic oscillations take 
place in the modified system 

Fig l,c demonstrates the dependence of the fourth 
largest Lyapunov exponents of coupled Rossler oscilla- 
tors 101 with the slight mistuning of the control param- 
eter uj d (uo d — 0.99) on the coupling strength e. Two of 
them, Af and A2 correspond to the behavior of the drive 
system, therefore they do not depend on e. Two other 
quantities A^ 2 are the conditional Lyapunov exponents. 
When the coupling parameter e is equal to zero, A^ > 
and A2 = 0. With parameter e increasing the second CLE 
A2 becomes negative (e rs 0.04), but the dynamics of the 
modified system © remains still chaotic (A™ > 0). With 
further increasing of e value the dynamics of the modified 
system © becomes periodical (see Fig. l,a,6), but GS is 
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FIG. 1: The bifurcation diagram (a) and the dependence 
of two Lyapunov exponent A™! (b) of the modified Rossler 
system @ on the parameter e. The third Lyapunov ex- 
ponent is about A™ ~ —9.7 and is not significant for our 
considering. The value of parameter e p when the modified 
system starts demonstrating the periodic dynamics is shown 
by an arrow. (c,d) The dependence of the Lyapunov expo- 
nent spectrum on the parameter e for slightly (ujd = 0.99) 
and greatly (cj^ = 1.3) detuned Rossler systems, respec- 
tively. The onset of GS is marked by an arrow. Con- 
ditional Lyapunov exponents are presented by dashed (X[) 
and dotted (A£) lines, (e) The bifurcation diagram for the 
non-autonomous modified system. The first equation in © 
is replaced by x m — —uj r y m — z m — ex m + Acos(Q,t) where 
A = 1.32, Q — 1.0 that simulates the parameters of the drive 
Rossler system. The value of parameter e p corresponding to 
the onset of the periodic oscillation is shown by an arrow 



not yet observed. It takes place only when the period- 
ical regime of the modified system © is stable enough 
(egs ~ 0.11). In this case the period-1 cycle is realized 
in the modified Rossler system. Below the critical value 
e c w 0.15 the modified Rossler system comes to the sta- 
tionary state. Note, that when the periodical regimes 
take place in the modified system, the value of the high- 



est CLE is slightly negative if GS is realized. As soon as 
the stationary state of the modified system © becomes 
stable the value of A™ starts to decrease rapidly. 

Note also, that the onset of GS is determined by the 
stability of the periodic regimes of the modified sys- 
tem @ which does not depend on the mismatch of 
the control parameters gd.r of the unidirectionally cou- 
pled oscillators. The stability of the periodical regimes 
is caused by the property of the modified system only. 
Therefore, the value of ecs should not depend greatly 
on the parameter mistuning (compare the values of egs 
for the uj d = 0.99 (Fig. l,c) and uj d = 1.3 (Fig. l,dVL This 
conclusion agrees well with numerical results of [18J. 

Let us briefly discuss why the onset of GS does not co- 
incide with any bifurcation point of the modified system 
(compare Fig. 1, a,b with Fig. 1, c,d). The cause of this 
non-coincidence is the influence of the external signal on 
the modified system. As it has already been discussed 
above the external signal (even if it is harmonic) can 
excite proper chaotic oscillations in the dynamical sys- 
tem with periodic dynamics. Therefore, the bifurcation 
points of the modified system under the external signal 
will be shifted in the direction of the large values of the 
e-parameter in comparison with the autonomous dynam- 
ics of the modified system. It is clear that the onset ecs 
of GS can not coincide with bifurcation point of the au- 
tonomous modified system. 

This statement is illustrated in Fig. 1, e where the bi- 
furcation diagram for the response system 10 under the 
external harmonic signal simulating the drive system sig- 
nal is shown. One can see that all bifurcation points of 
the modified system in the non-autonomous regime are 
shifted in the direction of the larger values of e-parameter 
(compare Fig. 1, a and Fig. 1, e). 

The same effects take place when GS is observed in 
the discrete maps. For example, GS takes place for the 
coupling parameter values e > Egs ~ 0-32 (see ^3| for 
detail) in two unidirectional coupled logistic maps 



Vn+i = f(y n ) + e(f(x n ) ~ f(y n )), 



(7) 



where f(x) — 4x(l — x). Following the concept described 
above one can construct the modified system 



z n+ i = (1 - e)f(z n ) = az n (l - z n ) 



(8) 



(where a = 4(1 — e)) and obtain that the value egs cor- 
responds to the value of a 2.72 of the equation ©. For 
such a value parameter an attractor of the logistic map 
is the stable fixed point x a — (a — l)/a. 

Let us briefly discuss now the case of GS between os- 
cillators of different types or when the coupling between 
oscillators is not diffusion. Several examples of such sys- 
tems are known (see, e.g., 0, Q). Obviously, if the 
coupling type is diffusion the difference of the system 
types does not matter and all reasons mentioned above 
remain true. But what happens when GS takes place in 
the systems coupled in different way, rather than in J3J)? 
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FIG. 2: Time realization y r (t) corresponding to the au- 
tonomous dynamics of the Lorenz system (solid line) and 
£GSUd(t) (dashed line) corresponding to the external signal 
introduced into the response system near the onset of GS 

One of the examples of such systems (see ^3[ for detail) 
is the coupled Rossler (drive) 

id = -a{yd + z d ), 

ijd = a(x d + ay d , (9) 
z d = a(p + z d {x d - c)) 

(a = 6, a = 0.2, p = 0.2, c = 5.7) and Lorenz (response) 

X r — &(llr *^r)i 

y r = rx r - y r - x r z r + ey d , (10) 

systems, where a = 10, r = 28, b = 8/3. It is known [l3[ . 
that the value of the coupling strength corresponding to 
the onset of GS is egs ~ 6.66. The amplitude of oscil- 
lations of the yd-coordinate of the Rossler system for se- 
lected parameter values is about 10, the amplitude of the 
^-coordinate of the Lorenz system being in autonomous 



regime (e = 0) is about 20. Obviously, the amplitude 
of the external signal ScsUd introduced into the response 
system (near the threshold of GS regime arising) is about 
60. So, the magnitude of the external force exceeds the 
level of proper system oscillations in several times. This 
situation is illustrated in Fig. 2 where the time series of 
y r (t) corresponding to the autonomous dynamics of the 
response system <|10(1 and the external force £GSUd(t) are 
shown. It is clear, that the great external force destroys 
completely proper dynamics of the response system, the 
phase trajectory of the Lorenz system is moved into the 
regions of the phase space with the strong dissipation and 
the mechanism discussed above causes the appearance of 
GS. 

In conclusion, we have explained GS appearance. The 
modified system approach has been proposed to demon- 
strate the reasons of GS arising. We have shown that 
the behavior of the response chaotic system is equal to 
the dynamics of the modified system (with the additional 
dissipation) under the external chaotic force. The cou- 
pling parameter increase is equivalent to the simultane- 
ous growth of the dissipation and the amplitude of the 
external signal. 
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